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SUMMARY 


The development of General Electric’s distortion methodology. Method D, 
has been documented. Method D was then used to correlate steady-state and 
unsteady data sets. A modification of Method D, known as Method D+ was 
documented and applied to a steady-state data set. 

A NASA-LeRC distortion methodology formulation, DIDENT, was investigated 
for accuracy of its correlation capability, ease of implementation, and appli- 
cability. Three different formulations of DIDENT were investigated. The first 
provided a correlation of loss of surge pressure ratio at constant corrected 
speed. The second formulation correlated loss of surge airflow at constant 
corrected speed. This was then combined with the first formulation to result 
in a correlation of loss of surge pressure ratio at constant corrected air- 
flow, The third formulation provided a direct correlation of loss of surge 
pressure ratio at constant corrected airflow. These three formulations were 
applied to steady-state data. 

One of the aspects of DIDENT is that a distortion pattern recognition 
technique is needed to properly choose the sensitivity parameter used in pre- 
dicting loss of surge pressure ratio or loss of surge airflow. Three different 
pattern recognition techniques were investigated. Two, developed during this 
work, are based upon radial and circumferential pressure defect terms. The 
third, developed by NASA-LeRC, is based upon maximization of the loss of 
surge pressure ratio. 

It was concluded that the formulation of DIDENT that provided a direct 
correlation of loss of surge pressure ratio at constant corrected airflow was 
the best DIDENT approach. This conclusion is based on applicability and use- 
fulness in all phases of engine design. This formulation of DIDENT provides 
a correlation quality that is comparable to Method D, 

Of the three pattern recognition techniques, the NASA-LeRC technique 
provided the best results. The first two techniques suffered because too many 
combined screen patterns were recognized as circumferential while the test 
engine saw them as radial patterns. 

The recommended formulation of DIDENT using the NASA-LeRC pattern tech- 
nique was then applied to an independent steady-state data set. The good 
correlation resulting from this indicates that DIDENT is a viable approach to 
distortion methodology. This formulation of DIDENT is not related to parallel- 
compressor theory. Because of this, to avoid confusion in the report, it will 
be referred to as Method E. 



An attempt was made to use the recommended formulation of Method E to 
correlate the unsteady data set used with Method D, Analysis of this work 
showed that the Method E sensitivities are functions of pressure defect level. 
In the general use of Method E, the levels of distortion used to define sensi- 
tivities must cover the range of anticipated distortion levels from flight 
tests. 



INTRODUCTION 


Engine inlet distortion has long been recognized as a source of aero- 
dynamic instability that can lead to engine surge. In recent years data from 
high performance inlet/engine systems have shown that engines can be sensitive 
to unsteady distortions that only exist for milliseconds. 

Various correlations have been developed to predict a given engine ' s 
tolerance to distortion. These have primarily been empirical in nature. 
Typically, pressure defect parameters are defined to describe the distortion. 
These parameters are then related to surge pressure ratio loss by defining 
various coefficients, or sensitivities. Using these sensitivities, then, for 
any given distortion pattern, a surge pressure ratio loss can be predicted 
to assess an engine’s stability with that pattern. Unsteady distortion analy- 
ses have typically been deterministic in that surge pressure ratio loss is 
computed as a function of time. It is assumed that the steady-state-derived 
sensitivities apply to the unsteady analysis. The one unknown parameter then 
is the data frequency content required to make calculated peaks in surge 
pressure ratio equivalent to known available surge pressure ratio. The indi- 
vidual unsteady pressures are filtered by analog or digital means to define 
the data frequency content. This process leads to the concept of critical 
time or critical frequency in that a time or frequency can be used to describe 
the filter characteristics. When the filter characteristics necessary to 
correlate data from a given type of engine or compressor are known, a set of 
unsteady data can be used to generate a prediction of surge pressure ratio 
loss. 


The foregoing approach to predicting surge pressure ratio losses is 
used with General Electric’s Method D, This methodology has been in use at 
General Electric since 1970. It has shown reasonable correlation capability 
for both steady-state and unsteady total pressure distortion data. This 
report describes the effort that generated Method D and presents examples of 
its application to both steady-state and unsteady data. 

A different concept of relating engine inlet distortion to surge pressure 
ratio loss has been developed by NASA-LeRC. This concept, called DIDENT, has 
its basis in the parallel-compressor modeling technique. The mathematical 
formulation is quite simple. Also, the coefficients, or sensitivities, 
required by the DIDENT formulation have the potential of being defined by a 
relatively small data set. The implementation of the DIDENT concept follows 
similar techniques used for Method D. 

Various formulations of DIDENT are analyzed, using steady-state data sets 
as a basis of verification. The best, most useful, formulation is then applied 
to a new data set for further verification. This formulation is also applied 
to a set of unsteady data. These results are used to discuss the applicability 
of the DIDENT formulation to engine inlet total pressure distortion analysis. 
Aspects of unsteady data analysis, such as analog or digital analyses and input 
or output filtering are addressed and the advantages of each are noted. 
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ANALYSIS APPROACH 


Some of the factors to be considered in the evaluation of a distortion 
methodology are: 

• The correlation capability of the methodology must be sufficiently 
accurate to allow confidence in its prediction. The current goal 

in the industry is an error band of +0,02 when the correlated parame- 
ter is surge pressure ratio loss. 

• The correlated parameter should be readily related to surge margin. 
Without such a relationship, it is difficult, if not impossible, to 
use a distortion methodology in stability assessments of engines 
early in their design phase. An accepted practice is to define loss 
of surge pressure ratio and surge margin at constant corrected air- 
flow. With these parameters, the designer is able to relate objec- 
tive operating lines, surge lines and distortion sensitivities prior 
to any detailed component design. This freedom is considered to be 
an important aspect of a distortion methodology’s usefulness. 

• The methodology should be easy to implement. That is, the definition 
of methodology sensitivities should not require an inordinately large 
data base. 

• The methodology should be applicable to any compression system (i.e,, 
fans, boosters or compressors). Further, it cannot show variations 
of sensitivity from engine to engine. For example, the sensitivities 
desired from a cell test must be applicable to flight tests using two 
different engines of the same design. 

With these criteria, then. Method D and DIDENT can be evaluated. Method D has 
been applied to steady-state data and to unsteady data. The unsteady data 
analysis includes a comparison of analog and digital analyses. A refinement 
of Method D is discussed and is applied to steady-state data. Three formula- 
tions of DIDENT are applied to steady-state data. The most applicable formula- 
tion is then applied to an Independent steady-state data set and to unsteady 
data. This unsteady analysis uses analog techniques. 

Five data sets have been used in this evaluation and are briefly described 
here. A more detailed description can be obtained from the references cited 
with each description. 

1. Data Set A - A set of 11 drift surges from a NASA-LeRC J85-13/axi~ 
symmetric inlet wind-tunnel test. The data are characterized by 
high turbulence levels, up to 8%, with no particular class of pattern 
(circumferential, radial) being predominant. These data are 
described in Reference 1, 


4 



2 . 


Data Set B - A set of 102 surges from a NASA-LeRC J85 -13 screen 
test. The data set consists primarily of circumferential, tip- 
radial, and hub-radial screens. Only 4 screens can be classified 
as combined. Reference 2 gives a description of these data. 

3. Data Set C - A set of 84 surges from a NASA-LeRC J85-13/screen test. 
The data set is characterized by a large number of combined patterns . 
Two engines were used within this set of data. Reference 3 gives 

a description of these data. 

Data Sets B and C are defined by a set of six equally— spaced rakes. 

Some screen positions were such that rakes were located at the edge of the 
screen rake. This tended to artificially suppress face average pressure and 
increase pressure defect extents. These rake values were corrected using 
NASA-LeRC-supplied correction factors during the DIDENT analyses. 

4. Data Set D - A set of screen test data from a variety of GE tests 
along with Data Set B, Data Set D was used in the development of the 
General Electric Method D, Reference 4 gives a description of these 
data. 

5. Data Set E — A set of 36 surges from a General Electric low press-ure 
compressor component screen test. The data are described in 
Appendix I. 

Two different definitions of loss of surge pressure ratio have been used. 
Figure 1 shows the relationships between the various terms for these two 
definitions. The loss of surge pressure ratio defined at constant corrected 
airflow (APRSW) is: 


APRSW = 


PRc - PRb 
PR c 


The loss of surge pressure ratio defined at constant corrected speed (APRSN) is: 


APRSN = 


pra - PRb 


A loss of surge airflow at constant corrected speed (AWS) is also defined as: 

w/e~ ^ w/e \ 

6 ^A 6 


AWS = 


w/e 


) 


5 


DISCUSSION OF RESULTS 


DEVELOPMENT OF METHOD D 


The General Electric Method D is an empirical approach to correlating 
pressure distortion to loss of surge pressure ratio. It was developed from 
other empirical methodologies (referred to here as Methods A, B, and C) . 

Method D and its predecessors are outlined in Reference 4. The equations 
used by Methods A, B and C and their implementation are included here to 
provide a background for the development of Method D. The methodologies A, 

B, C, and D are defined in Figures 2 through 5. Method A, shown in Figure 2, 
uses a Fourier series to describe each ring’s circumferential pressure varia- 
tion. Each ring is weighted inversely to diameter so that hub circumfer- 
ential distortion contributes more heavily than tip circumferential distor- 
tion. Radial distortion is a function of ring average pressure, again 
weighted inversely to diameter. Pressures are normalized by inlet dynamic 
pressure, q, in an attempt to collapse speed effects. Complex patterns use 
a direct linear superposition of circumferential and radial contributions. 

Method B is shown in Figure 3. All pressures are normalized by face- 
average pressure. Here the circumferential variation in pressure is described 
by a Fourier series. The series coefficients are summed and ring-weighted 
to define a circumferential index, A(^. This is made equivalent to surge 
pressure ratio loss, APRS, for pure circumferential test data: 


^C AVAIL 


Thus, the circumferential parameter 

™c = 

is equivalent to APRS.„.,^ . 

AVAIL 

The radial distortion parameter is similarly defined. Ring average 
pressures are described by a Fourier series. The series coefficients are 
summed and weighted to define a radial index , Ar. This is related to APRS 
for pure radial test data: 

and 

= 'Si * \ 
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Complex patterns are treated with a linear superposition factor of 1.0: 

ND^ = ND^ + NDj^ 

Note that ND^ is equivalent to APRS AVAIL* 

Method C is shown in Figure 4. This technique relies on level param- 
eters rather than Fourier coefficients. The circumferential ring level 
parameter is: 


IDC 


P . - P . 

ring avg ring min 


P 


Face avg 


This level parameter is modified by shape, extent and multi-per-rev function 
(1/n)^, factors to define a sensitivity. This is, in turn, summed and 
weighted to develop a Radial data are similarly processed. The 

radial level parameter is defined: 


IDR = 


Pt. - P ^ 

Face avg ring avg 
p 

Face aVg 


This parameter is weighted with extent factors k and separate tip and hub 
sensitivities are defined. Complex patterns are processed by using a linear 
superposition factor of 1,0. The final parameter, APRS^otal equivalent 
to aprs^y^il. 

Method D, which is an extension of Method C in that pressure defect 
definitions are similar, is shown in Figure 5. Here APRScircum for each 
ring is defined by the same level parameter used in Method C, modified by a 
"multi-per-rev" factor, a shape factor, an extent factor and a sensitivity. 
APRScircum then chosen as the maximum of the average of the two tip ring 
APRS values or the two hub ring APRS values. 


Radial data are handled by an extent coefficient and a sensitivity. 
They are, in general, different for the hub and tip regions. The largest 
APRS is chosen as representative of stall margin loss due to radial distor- 
tion, Method D is unique in that complex patterns are processed using a 
nonlinear superposition factor. It is a function of the ratio APRSjjad/ 
APRS^,^ and operates on the APRScircum term: 


APRS^ , = b 
Comb 


* APRS , + 

circum 


APRS 


RAD 


Again, APRS^oiab fs equivalent to APRS avail a prediction of surge 

pressure ratio loss is computed. 
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It is estimated in Reference 4 that Method D requires 22 patterns at 
each speed to conpletely define the sensitivities. 

The correlation of each methodology with the General Electric Data Set 
D, which includes Data Set B, is shown in Figures 6 through 14. Only the 

bounds of Data Set B are shown. Note that Method A derives a correlation 

with 4PRS^y^jLj that is related only to a given data set. The parameter 
is not universally related to APRSjs^VAIL* this reason, Method A data 

are shown on individual figures for each subset of data in Data Set D. This 
is not the case with Methods B, C and D. All of these derive a number 
(NDT or APRSc^2.c^ that is universally related to APRS^^y/^jL, regardless of 
the data set used. For these methodologies, the entire Data Set D is 
presented on one figure. 

A qualitative cou^arison of these correlations can be made through the 

use of the standard deviation of each subset. This is shown in Table I, 

derived from data in Reference 4. The standard deviation, a, is defined as: 


N .2 

1=1 

N - 1 

where N is the number of points 

X is the methodology parameters NDT or APRSq^-j^^ . For K^, the 
standard deviation is that obtained from a straight line fit 
of vs APRSxest* 


Table I. Standard Deviation. 


Subset 

A 

B 

C 

D No. 

of Points 

GE4 (72T, 85T) 

0.0345 

0.0244 

0.0141 

0.0073 

9 

GE4 (72P, W85, Va5) 

0.0138 

0.0168 

0.0237 

0.0167 

13 

GE-FIOO BU/6 

0.0295 

0.0342 

0.058 

0.0190 

19 

GE-FIOO BU/7 + IGV 

0.012 

0.0197 

0.0075 

0.0067 

13 

GE-FlOO BU/7 + IGV + OS 

0.020 

0.0052 

0.0074 

0.011 

18 

NASA-J85 (set B) 

0.04 

0.0268 

0.0234 

0.0185 

100 

Overall 

0.0337 

0.0248 

0.0270 

0.0165 

172 

data show that Method 

D provides the best correlation on 4 

of the 6 


subsets and is the best method of those shown on the overall basis. 

The development of the empirical constants used in Method D for the 
NASA Data Set B are shown in Figures 15 through 23. Circimiferential 
sensitivies Kq, are the slopes of the straight lines in Figure 15. These 
parameters define Kq as a function of corrected speed and extent of the low 
pressure region. These values are tabulated in Table II. 
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Table II. 

Data 

Set B Values of Kq. 


Angular Extent, 

Degrees 

Corrected Speed 

215 

141 

100 

87% 

2.32 

1.65 

1.10 

90% 

3.18 

2.39 

1.90 

93% 

2.40 

1.50 

1.37 

97% 

1.00 

1.34 

1.88 

100% 

1.4 

0.93 

0.93 

Screen Angular Extent 

180° 

0 

o 

0 

c 


The above values were plotted as a function of extent in Figure 16. The 
value of Kq at 20° extent is assumed to be zero. These curves are shown 
as straight line segments between data points in accordance with the General 
Electric Distortion Analysis Program (GE DAP) method of interpolation between 
data points. Note also that the value of K(^ at 180° has been interpolated 
using faired curves, rather than straight line segments between 140° and 
216°. 


The extent function, EF, is derived from these data since 

Kc) 0 

K„) 180° 

The extent function is plotted in Figure 17. These values of extent function 
were input to the GE DAP program along with at 180° to completely define 
the circumferential sensitivity. 

The "multi- per- rev" function, MR, was defined from data shown in Figure 
18. The data are limited and the scatter is large so that the "multi-per-rev" 
function was chosen to be 1.0 for all speeds. 

Radial sensitivities are defined by two parameters, Kr and a. These 
are interrelated through the equations : 


APRSj^p = KRjjP [AP/P)5 + * AP/P)4] 

~ [AP/P)x + “[jjjj * AP/P)23 

APRSp^ = max of (APRSjpp, APRSjjjjp) 
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^RTIP* ^RHUB ^HUB were computed from the first two equations using 

a linear least squares multiple regression analysis with APRS)test _ APRS)|-a(j. 
The resulting values of these four coefficients were used to construct the 
lines in Figures 19 and 20, The data in these figures show a reasonable fit 
with the constructed lines. 

The data from these figures were then used to define the variation of 
Krh, Krj., ajj and with speed. The results of this are shown in Figures 

21 and 22. 

The preceeding circumferential and radial sensitivities were used 
to define APRSciRC ^nd APRSj^j^) for the various combined patterns in the 
data set. These were, in turn, used to develop the superposition factor, 
b. b is defined as: 


APRS 


b = 


AVAIL 


- APRS 


RAD 


APRS 


CIRC 


It is a function of 


APRS 


RAD 


APRS 


CIRC 


and corrected speed. 


The data shown in Figure 23 define the b-f actor for these data. Note 
that they are not a function of speed in this case. The addition of any radial 
distortion makes the J85-13 compressor less sensitive to the imposed circum- 
ferential distortion. 


The correlation of all points in Data Set B, using these constants, is 
shown in Figure 24. From Table I, the standard deviation is 0.0185. 


REFINEMENT OF METHOD D 

In 1972, General Electric evaluated Method D with the objective of 
refining its correlation capability. This evaluation was based upon the 
Data Set C. The result of this work has become known as Method D+. As with 
Method D and its predecessors, Method D+ is an empirical approach. The 
primary difference between Methods D and Df is the inclusion of radial 
weighting terms in Method IH-. The formulation of Method D+ is as follows: 

• Circumferential 

P . - P . 

ring avg , ring min . 

i is ith ring 

face avg 
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1 ^ AP\ . ^ ^ AP\ 

1 J ^ ^ ^ rings where ^ ) gi 

C " i ^ "Ci ^ '^i 


AP/P^ 


= X CEP X REF X MRP 


E V.P y 




^^^CIRC = *^c* 


(f) 


Radial 


P. - P . 

face av rxng av^ 




face av 


i is ith ring 


AP/P) :>. 0 


300° or more of ring must be below P 
for AP/P^^^ to be > 0. 


face av 


^ “i * 


Ri 


"i AP/P 


) 


Combined 




^ f (^P^BAD> 




The application of Method VH- to Data Sets B and C has been limited to 
87, 93, and 100% corrected speeds where the bulk of the test data are located. 

The circumferential data are analyzed by first plotting (AP/P)q versus 
Kc, where = APRS^Y^l/ (AP/P)^ for 180®, full-span screens. This develops 
the basic circumferential sensitivity, to which all circumferential data 
are normalized. The first normalization considers angular extent. Full span 
screens of extents less than 180* are plotted as the 180* screens were. This 
is shown in Figure 25 for the three speeds analyzed. The ratio of sloues 
between 180* and any other extent is the circumferential extent factor, CEF, 
which is plotted versus angular extent in Figure 26. The second normaliza- 
tion considers radial extent of circumferential screens (the partial extent 
tip and hub radials in Data Set C, for example). Here, Kq is plotted versus 
C^/P)C X CEF. 

The result, shown in Figure 27, allows a definition of radial extent 
factor, REF, as the ratio of slopes between the radial extent = 1 and any 
other radial extent. REF is plotted versus radial extent in Figure 28. This 
defines the circumferential screen sensitivities, all normalized to 180*, 
full span sensitivity. 

Radial data are considered on a ring-by-ring basis. That is, a sensi- 
tivity, ai = APRSAVAIL/(^P/I’)Ri plotted versus (AP/P)Rj;^ for each ring. 

Data where more than one ring has a contribution [(AP/P)r£> 0] are considered 
by assuming that the rings contribute to APRS linearly. For example: 


where i is ith ring - the only ring to contribute to APRS in this data point 
k is kth data point 

where i is same ring as above 
j is another ring 
k is same data point as above 

1 is new data point where ith and jjA rings contribute to APRS. 

This process can be used to completely define the ct^’s as shown in 
Figure 29. The gain in APRS with hub radial screens causes some negative 
values of . 
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The circumferential and radial data so far analyzed can now be used 
for combined patterns, using a superposition function: 


A = 




where 


^^^CALC ■" ^ ^^^^RAD 


The superposition function was determined by solving the above equation 
for A and plotting these values versus APRSj^, as shown in Figure 30, 

Note that, in order to calculate APRS for circumferential-only data, 

A " 1 at APRSg^ = 0. 

There is a high degree of scatter in the data in Figure 30. The curves 
drawn through these data are considered the best fit to the data. 

Putting all of these sensitivities together produces the correlation 
shown in Figure 31. Here APRSq^q is plotted versus APRS^y^jL . There are 
some rather large deviations from the main correlation band. While they are 
limited in number, the underlying data must be considered valid and these 
points represent a failing of Jfethod D+. 

The standard deviation is 0.025 for this work. However, when comparing 
Methods D and D+, the comparison must be on an individual data set basis 
because Method D coefficients were developed on the individual data sets. 
This is shown in Table III where the various standard deviations are tabu- 
lated. 


Table III 


Comparison of Methods D and D+ 


Method D 


Method IH- 


Data Set B 0.018 
Data Set C 0.036 
Combined 0.027 


0.011 

0.039 

0.025 


This table shows that Method D4- provides a better correlation for the Data 
Set B, largely pure circumferential and radial patterns. However, when 
combined patterns are introduced with Data Set C, it is no better than 
Method D. 

The conclusion has been reached that Method D+ is not preferred over 
Method D because it provides a correlation no better than Method D and 
requires more detailed information (i.e., more screen patterns) than Method D. 


I 
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APPLICATION OF METHOD D TO UNSTEADY DATA 


This section discusses the analysis of NASA-LeRC Data Set A using 
Method D as programmed in the General Electric Djmamic Distortion Analysis 
Program (DAP). The analysis technique used in DDAP processes the unsteady 
pressure data time^'Slice by time-slice to compute surge pressure ratio loss 
(APRS) via the Method D equations. The output, then, is APRS as a function 
of time (DDAP can plot up to six parameters on up to six plots for a total 
of 36 parameters versus time) . 

Input pressure data are digitally low^pass filtered using a sliding 
window average. The aiqplitude characteristics of this type of filter are 
shown in Figure 32 where this filter is compared to a 5— pole linear phase analog 
filter. They are seen to be quite similar in character. Note that two 
characteristic frequencies are defined in Figure 32 - the 3 dB down point, 
and the pole, or zero aiq)litude ratio point. The pole frequency is not 
realizable in an analog filter and conq>arisons between analog and digital 
filters are best made using the -3 dB frequency. However, in DDAP > ^nd 
in this report, the pole frequency is the inqilied frequency of the filter 
in that : 


Averaging Time = — ^ 

pole 

For the filter employed in DDAP the relationship between the pole and -3 dB 
frequencies is: 

^(-3 dB) "0*^^ ^(pole) 

This input filtering is the control over the time dimension of the unsteady 
analysis. Because all Method D sensitivities and coefficients in DDAP are 
just those derived from steady-state screen data (Hata Set B) , the input 
filter is the only Independent variable that controls the amplitude and 
frequency characteristics of APRS, This leads to the concept of "critical 
time," Each data point, consisting of a digital record of pressure time 
histories leading to engine surge, is successively processed with different 
filter averaging times. The APRS time histories for each averaging time 
will typically show a peak just prior to surge. The anqilitude of this peak 
is compared to the known surge pressure ratio margin (^PRS^Y^l) . The 
"critical time" is just that filter averaging time that causes the peak 
APRScalC prior to surge to be equal to APRS^V^AXL* When a group of data 
points are available for analysis, the "critical time" is chosen so that 
the standard deviation of APRSq^q compared to APRS^y^jL minimized. 

The DDAP analysis has been used with Data Set A and Method D sensi- 
tivities derived from steady-state analysis of Data Set B. A t 3 rpical time 
history of loss in surge pressure ratio (APRScALC) is shown in Figure 33. 

The measured value of APRS^^y^^L peak in APRScALC that induced surge 

are shown. This figure shows an apparent shift of the data in time as a 
function of averaging time. This is most easily seen by examination of the 
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peak in APRSq^q that occurs after compressor surge. This apparent shift 
is due to the way in which the averaging is implemented. The value of any 
individual pressure (and, therefore, APRS^^q) at a given point in time is 
the average of that pressure from the given point in time to the leading 
edge of the average window. In equation form; 

tj_ + N-1 

I S 

t . 

1 

ti is time slice being analyzed 
N is time slices being averaged 
Pj is individual values of pressure at each 
time slice. 

In general, there would be an additional phase shift due to the number 
of time slices advanced between averaging processes. For all data reported 
herein, the advance was one time-slice. It has no effect in these data. 

This apparent shift in time has no effect on the analyses performed on 
Data Set A, but in comparisons of wave shapes, it must be considered. 

In Data Set A, the location of surge was chosen as that point where 
face average began to rise above the previous fluctuations of face average. 
This is the first indication of the surge overpressure. An analysis window 
was then defined from this surge location. The highest peak of APRS within 
this window represents the peak stall pressure ratio loss for a given 
averaging time. The peak, in general, moves within the window as averaging 
time is varied. That is, one unique peak in APRS is not maximized for all 
averaging times. This is because the wave shape of APRS may' be, in part, 
composed of nearly simultaneous minima from different probes. Input fitering 
then, will affect the waveshape of APRS through the filtering effect on 
individual probes. In any analysis of unsteady distortion, the criteria used 
to define the peak APRSc^C niust be defined. The approach used for these 
data is to assume that the engine surged due to some distortion that created 
a APRS of sufficient magnitude and duration to induce the surge. The time 
required to propagate a disturbance from the measurement plane to the com- 
pressor discharge, allow one compressor revolution, and propagate a distur- 
bance forward to the measurement plane is conservatively estimated to be no 
greater than 20 milliseconds. The peak APRS that induced surge, then, 
should occur no further than 20 milliseconds prior to surge. Peaks in APRS 
that occur at times greater than 20 milliseconds prior to surge have no 
effect on engine surge, regardless of their amplitude or duration. They 
represent an error in the correlation capability of the methodology used in 
the unsteady analysis. The peak values of APRS, obtained within 20 milli- 
seconds of surge, are shown in Table IV. Averaging times are normalized by 
the time required for one coiq>ressor rotor revolution at its design speed 
of 16,500 rpm. 
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Table IV 


Data Set A Surge Pressure Ratio Losses 
Point AVAIL -^^^REV «25j^y APRS) ^ 3 ^^ APRS)i^y APRS> 2 j^y 


82 

0.147 

0.168 

0.15 

0.145 

0.137 

0.132 

103 

0.156 

0.248 

0.202 

0.215 

0.162 

0.055 

141 

0.079 

0.16 

0.13 

0.09 

0.055 

0.027 

148 

0.183 

0.293 

0.25 

0.18 

0.111 

0.078 

154 

0.216 

0.303 

0.237 

0.15 

0.113 

0.065 

162 

0.238 

0.355 

0.332 

0.294 

0.263 

0.241 

164 

0.237 

0.166 

0.160 

0.153 

0.141 

0.140 

200 

0. 066 

0.147 

0.118 

0.078 

0.050 

0.043 

214 

0.009 

0.113 

0.107 

0.103 

0.087 

0.073 

219 

0.008 

0.132 

0.128 

0.118 

0.097 

0.057 

261 

0.244 

0.282 

0.268 

0.233 

0.196 

0.177 


The data shown in Table IV were used to develop a histogram of the 
number of surges where: 


APRS 


CALC) 


^ APRS 


A.T. 


AVAIL 


APRS 


CALC) 


A.T. + AA.T. 


This is shown in Figure 34. The peak of the histogram is at 0.25 to 0.5 
rotor revolutions, but it is not conclusive that this time is the critical 
time of the J85--13 compressor. The histogram shows data points that do 
not agree with the basic distribution of the histogram. Referring back to 
Table IV, these data can be grouped and qualified as follows: 


o Points 82, 103, 141, 154, 200 and 261 are described by APRS). 125 > 
APRS avail ^ APRS) 2 . These points indicate that the critical 
time is between 0.125 and 2 rotor revolutions. 

o Point 164 is described by APRSaVAIL > APRS).i25* Because of the 

eight points mentioned above, and the low senitivity to 30" extent 
per rev screens, it is highly unlikely that a realistic critical 
time less than 0.125 rotor revolution exists within these data. 
This point is interpreted as one containing an error in either 
APRSavAIL or APRScalC» or both. It is not considered further in 
this analysis. 

o Points 162, 214, and 219 are described by APRSavAIL APRS) 2 . 
Because of the extremely low APRSavAIL ^or 214 and 219, either 
the engine should have surged using a steady-state analysis or the 
APRSavAIL error. These 2 points (214 and 219) out of the 

11 drift surges from Data Set A will not be considered further. 
Point 162 would indicate a critical time slightly greater than 2 
rotor revolutions. This point is eliminated because of the 
preponderance of the data lying well below 2 rotor revolutions. 
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With these "outlayers" eliminated from the data, the critical time 
can be investigated. This has been done by computing the standard deviation 
of the data set at each averaging time analyzed. The standard deviation 
compares the peak to ^PRS^VAIL each averaging final. 

The results are plotted in Figure 35 as a function of rotor revolutions. 
The data show a definite minimum in a at 0.5 revolution. Because points 
162 and 164 cannot be eliminated from the analysis on grounds other than 
they do not show a critical time between 0.125 and 2 revolutions, the 
standard deviations were computed including these points. The results are 
also shown in Figure 35. Again, the standard deviation shows a minimum at 
0.5 revolution. 

The value of the standard deviation at 0.5 revolution is a measure of 
the accuracy of the methodology used in correlating the unsteady data. From 
Figure 35, a = 0.037 at 0.5 revolution, while a - 0.0185 for the steady-state 
Data Set B. This would indicate an accuracy for Data Set A of ~0.5 that for 
Data Set B. However, the difference in the nimiber of points within each 
data set could cause these standard deviation comparisons to be misleading. 
Figure 36 shows the basic correlation of with APRSqalC Data 

Sets A and B, The data from Data Set A are seen to correlate quite well 
except for two points. It can only be concluded that Data Set A is not 
large enough to statistically compare its fit to Data Set B. Figure 36 
would indicate that the quality of the fit of Set A is equivalent to the 
fit obtained with Data Set B. 

Figure 36 also demonstrated an important difference between Data Sets A 
and B. Values of APRS for Data Set A are sometimes larger than those for 
Data Set B. This is important because Set B defines the sensitivity coeffi- 
cients used in DDAP to analyze Set A. There is no assurance that these 
sensitivities may be accurately extrapolated to higher values of distortion 
as has been done in this analysis. However, referring back to Table IV, it 
can be seen that 0.5 rotor revolution will still provide the minimum standard 
deviation for those data points in Set A with APRS^Y^jL values within the 
bounds of Data Set B. 

It can be stated that the analysis of Data Set A shows that the J85-13 
has a "critical time" of 0.5 rotor revolution using the following assumptions 

• The sensitivities defined by Data Set B can be extrapolated to 
the higher distortion levels of Data Set A. 

• The peak APRSqaLC that induced surge occurred within 20 milli- 
seconds prior to the surge. 

• The concept of "critical time" is valid. That is, sensitivities 
based upon steady-state data can be applied to unsteady data. 


ANALOG VS. DIGITAL ANALYSIS 


The analog analysis of Data Set A was performed in General Klectric’s 
Stability Measurements Analysis Laboratory (SMA.L) . A photograph of the 
equipment used for the analysis of Data Set A is shown in Figure 37. The 
computers were programmed as shown in Figure 38. The objective of the 
analysis at the time it was done was to screen the data records for peaks 
in the Method D parameters AP/P)^ and AP/P)jj^. This analysis will form the 
basis of the discussions in this section. 

Referring back to Figure 38, the analog data were low-pass filtered at 
the playback discriminators to 200 Hz (-3 dB) and then summed with steady- 
state values to make up the 30 filtered time-unsteady pressures. These were 
then processed to compute the various parameters shown in Figure 38. 

Figure 39 shows typical digital and analog analyses for a point in Data 
Set A. The data were processed with comparable input filters. Each figure 
contains 250 millisecond records of AP/P)Cmax» AP/P)x ^^y and AP/P)cmax 
AP/P)Ruiax* Visual inspection of these wave forms show that the wave shapes of 
digital and analog analyses are quite similar. One measure of the similarity 
is a comparison of peak values. This is shown in Figure 40 where a peak value 
for all three of the above parameters in each of nine digital records is 
plotted against the corresponding peak value in its companion analog record. 
The difference between analog and digital analyses is seen to be ± 0.01 of 
AP/P)c or AP/P)r. 

When comparing cost and time required for analysis, some assumptions 
must be made as to the starting point. It will be assumed that the digital 
computer use is charged to the analysis, the analog computer use is free, 
the data are on an analog magnetic tape, and that 10 data points are to be 
processed to determine peak APRSGjy_,c to stall. With this starting 

point, the following table of required steps and the estimated cost of each 
may be made. A rate of $25/man hour is assumed. 


Table V 


Comparison of Unsteady Data Processing Costs 


Digital j 

1 Analog 

Item 

Manhours 

Cost 

Item 

Manhours 

Cost 

Digitize Data 

20 

500 

Prog. Analog 

16 

400 

Qualify Digital 
Data 

10 

*450 

Checkout 

8 

200 

Process Data 

20 

*1200 

Process Data 

20 

500 


Total 

$2150 


Total 

$1100 


* Includes digital computer usage charges 

* Based on assumed rates and times 
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This table serves to estimate costs. Overall timing must be estimated 
separately due to turnaround time in the digital computer. The estimate 
is two points per day for qualification and processing of digital data, run 
in parallel. Digitation could be performed in two days, qualification and 
processing would require six days, for a total of eight days to process the 
data. Analog analysis is a serial procedure through the steps in Table V - 
a total of 5-1/2 days. Analog analysis is about 1.4 times faster and 2 times 
cheaper than digital analysis. Assuming that digital data are "truth", analog 
analysis has an error of ±0.01, based upon the analysis of Data Set A. 

With all data analysis efforts, engineering judgement is necessary to 
determine the best procedure. However, a good philosophy is to perform as 
much analysis as possible on the analog system and resort to digital analysis 
only when its capabilities offset its increased cost. 


INPUT VS. OUTPUT FILTERING 


The effect of input as opposed to output filtering on APRS can be seen 
by considering the Method D equations as programmed in the Dynamic Distortion 
Analysis Program (DDAP) . 
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If the same equation is output filtered the form is: 
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(Eqn. 2) 


p , = p . + p . 

«iin min min. 
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where subscript ss is steady-state component and (') is unsteady component. 
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(Eqn. 3) 


where E of first order P 


terms are zero by stationarity. 


Equations (l) and (3) are different in formulation because summations are not 
distributive. In Equation ( l) , individual terms are averaged over N samples. 

In Equation (2) N ratios are averaged. While this certainly effects the 
results, perhaps the largest difference in input versus output filtering is due 
to the treatment of minima. Equation (l) uses a minimum from one probe whose 
average over N samples is minimum. Equation (2) uses a minimum that is the 
average of N minima - not necessarily from one probe. One would expect that 
the average of N minima will be less than the minimum probe averaged over N 
samples. Ihat is: 



This inequality is demonstrated in Figure 41 where AP/P)^ for each ring Is 
plotted for data point 148. Output averaged data are seen to be equal to or 
greater than input averaged data. 

A similar relationship may be obtained for the Method D radial parameter > 
AP/P)r. That is: 
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However, because of the nonlinear relation of the basic circumferential 
and radial distortion parameters with the inequalities noted above 

do not hold when comparing APRScALC input and output filtering. 

This is shown in Figure 42 where ApRSc^LC plotted for data point 103, 
Here APRSc^c)q is generally, but not always, greater than APRSc^c)^;, The 
only conclusion to be drawn is that input and output filtering are not equivr^ 
alent. Further, the proper technique is input filtering because the basic 
premise is that the engine is responsive to pressure defects that exist for 
some length of time. Input filtering is the control over that length of time. 


DIDENT FORMULATION 


One of the disadvantages of most of the distortion methodologies currently 
used is that a significant number of screen patterns must be tested to define 
the methodology coefficients. An approach developed by NASA-'Lewis Research 
Center (LeRC) called DIDENT (Distortion Identity) shows promise in that, 
potentially, the number of screen patterns can be significantly reduced and 
still correlate complex patterns. A description of the approach is given in 
Reference 5. 

DIDENT is based on the parallel-compressor model and, for the J85-13 
data analyzed herein, assumes that the compressor discharge total pressure 
is uniform for any inlet distortion pattern. Following the parallel- 
compressor model then, equations have been developed in Reference 5 that 
relate distortion level to loss in surge pressure ratio at constant corrected 
speed (APRSN) . These equations are; 


APRSN 


AVAIL 


1 - 


(P,/P,)D 


N / = cons t , 


K. = 
1 


H ^3/^^min,60°^2,r ^ D 


(P3/P2) 


N/y^ = cons t . 


APRSN 


CALC 


(Pmin,60®)9 
= 1 - ^ 


X K. 
1 


P 3 - discharge total pressure 
P 2 “ inlet average pressure 


D = distorted surge line 
C =* clean surge line 


(P - - inlet minimum pressure averaged over 60 * 

min, 60 ® 2 r- • j 1 c 
in rings 1 or 5 

Kj[^ - sensitivity, a function of speed, pattern and engine. 
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Reference 5 has indicated that pattern effects are limited to circum- 
ferential, hub- or tip- radial effects. Thus, three separate curves, Kc, Kg, 
and K;r can be defined as speed functions for each engine. These curves, 
from Reference 5, are shown as Figure 43. Note that = 1.0, just the value 
predicted by parallel- compressor theory. 

The data in this curve show a strong engine- to- engine variation of the 
hub- and tip-radial sensitivities. This effect of an individual engine on the 
DIDENT sensitivities will be noted for all correlations of the DIDENT approach 
in this effort. It is discussed later in the report. 

A problem in using the DIDENT approach is the development of a pattern 
recognition technique to understand which Kj[ should be employed for any 
given complex pattern. 

In order to compare various pattern recognition techniques, it is 
necessary to know the data scatter inherent in DIDENT resulting from scatter 
in the Ki curves. This has been done for the APRSN correlation by manually 
choosing the Ki for each pattern. The APRS scatter is due only to scatter, 
since we have an identity. 

The results are shown in Figure 44 for Data Set B and Data Set C. The 
standard deviations of each correlation are noted in the figure. With this 
figure as a baseline correlation, pattern recognition techniques may be 
defined and investigated as to accuracy. 

Two different techniques have been investigated in this analysis. The 
first is based upon the radial pressure defect terms in Method D. This 
technique is: 




> 0 




Face 


> 0 




P 

Face 


The maximum value of AP/P)q, AP/P)rji, or AP/P)g selects the corresponding 
value of K^, Kj', or Kg. Data Sets B and C, correlated with this technique, 
are shown in Figure 45. The Data Set B, with its many pure patterns, is 
reasonably correlated. Notably, the two screens that differ from the base- 
line correlations are both combined patterns. Screen 24, a partial- extent 
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tip- radial pattern, was seen as a circumferential pattern at all speeds by 
this pattern recognition technique. The engine reacted to screen 24 as If 
It were a tip- radial pattern. Screen 26, a tip radial plus one-per-rev 
circumferential, was viewed by the pattern recognition technique as a tlp- 
radlal. The engine reacted to It as a circumferential screen pattern. 

Investigation of Data Set C shows similar results. Screens 17, 18, 19, 
21, 23, and 25 were viewed by the pattern recognizer as circumferential 
patterns. The engine reacted to screens 17, 18, 19 and 21 as hub-radlal 
patterns - they are all partial hub-radlals. The engine reacted to screens 
23 and 25 as tip- radial patterns - they are both partial tlp-radlals . 

The first technique, then. Is not capable of resolving all of the 
patterns In Data Sets B and C such that the loss of surge pressure ratio at 
constant speed Is predicted for all patterns. 

The second technique employed to discriminate the pattern character- 
istics was designed to view combined patterns as If they were pure circumfer- 
ential or radial patterns. This technique uses the following equations: 


^radj^ Pmln» 

Face 

^radj^ 

— 0 1 Is 1 th ring 

Face 

Pradi = 2 PRingAvg^ ' 50^ 

^min* 60 minimum pressure integrated over 60® in i th ring 

1 

PRingAv^ is ring average pressure of ith ring 
Ppace face average pressure 

With these definitions, AP/P)q and AP/P)j^ reduce to those used in Method D 
for pure circumferential and radial patterns, respectively. Combined patterns 
have either AP/P)c or AP/P)r suppressed, depending upon details of the pattern. 
This tends to make most combined patterns appear "pure." 

These equations were employed to compute the above circumferential and 
radial parameters in each ring. The circumferential parameters were then 
summed as follows: 
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was then used to select Kq, Kjj, or Kip, 


The resulting correlations for Data Sets B and C are shown in Figure 46 . 
Technique number two differed from the baseline correlation similarly to 
technique number one. The same screens noted for technique number one were 
viewed in the same way as technique nun±>er one did. In addition, technique 
number two viewed screen 23, Data Set B, as a circumferential pattern while 
the engine reacted as if it were a hub-radial pattern. Screen 23 is a 
partial extent hub-radial pattern. It also saw screens 15 and 24, Data Set 
C, as circumferential patterns while the engine reacted to them as if they 
were tip-radials. Both screens are partial tip-radials, with screen 15 
having two-per- revolution 90° extent tip sectors and 20° extent circumfer- 
ential sectors; screen 24 is a 120° one-per- revolution tip sector. 

The second technique proves to be less acceptable than the first. This 
technique did as anticipated - it viewed many combined screens as a pure 
pattern. However, too many patterns were viewed as pure circumferential 
rather than pure radial. This could be modified by changing the definition 
of AP/P(^. This was not done because, as later discussed, the utility of 
APRSN as an engine design program is limited. 

In general, the engine reacts to partial radial patterns as if they were 
pure radial patterns. The engine reacts to combined patterns with full-span 
circumferential components as if they were circumferential. Neither tech- 
nique was capable of fully resolving the patterns in Data Sets B and C in this 
manner. Of the two techniques investigated, the first is best for APRSN 
correlation. 


CORRELATION OF LOSS OF SURGE AIRFLOW 


As a continuation of the work on correlating Data Sets B and C using loss 
of surge pressure ratio at constant speed, a correlation of loss of surge 
airflow at constant speed has been obtained. These two correlations were 
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combined to predict the loss of surge pressure ratio at constant corrected 
airflow. Both pattern recognition techniques previously discussed were 
carried through these correlations . 

The method used in correlating loss of surge airflow is very simple, 
but was arrived at only after many other methods were investigated. These 
will be discussed, but no complete correlations of these attempts have been 
presented. 

The first approach was to utilize the parallel- compressor model as used 
in the basic DIDENT formulation. The sketch in Figure 47 for a two-sector 
parallel compressor shows the typical results of this effort. Here, the 
predicted airflow is that point on the straight line that crosses the pre- 
dicted pressure ratio at surge. The calculated surge airflow was generally 
higher than the actual surge airflow. Further, this calculation would 
become rather complex for multisector parallel compressors required by 
complex patterns. The approach was discarded due to inaccuracy and complexity. 

The second approach is based upon an observation that distorted speed 
lines may be derived from clean speed lines through the following technique. 
Define the coefficients: 



R 


2 



ri^ = Max 


Plot and R2 versus distortion with corrected speed and pattern as bivariant 
parameters and fit curves through these data. These curves then give values 
of Rl and R2 for given distortion, pattern and corrected speed. These can 
be used to compute a distorted airflow and pressure ratio based upon the 
known clean surge airflow and pressure ratio. The computed values do tend to 
lie on the measured distorted speed line, but there was no correlation with 
the distorted surge airflow and pressure ratio. This approach was discarded 
as too inaccurate. 

At this point, strictly empirical approaches were tried. They included 
curve-fits of the distorted maps which were too con5)lex and inaccurate. The 
auproach finally used, which resulted from the curve-fit work, was to simply 
define a loss of surge airflow as: 
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and plot it versus 1 - (Pm-f n » 60) 2> r/f2 constant corrected speed and 
pattern type (circumferential, hub-radial, tip- radial). These plots are 
shown in Figures 48 through 50. 

These correlations show that the J85-13 cojiqjressor tends to lose airflow 
with hiib-radial patterns and gain airflow with circumferential, tip-radial, 
partial- hub -radial, and partial tip-radial patterns. Further, Data Sets B 
and C give different correlations, indicating a difference in engine- to- 
engine behavior which was also noted in the DIDENT sensitivities of Reference 
5. 


The curves shown in Figures 48 to 50 were then used to correlate 
AWSavAIL with AWSq^q. This was done by entering each curve with known 
values of distortion level, pattern type (as determined by identif icatj on 
techniques numbers 1 and 2, discussed earlier) and corrected speed. These 
allowed a prediction, AWScALC to be determined. The resulting correlations 
of AWS are shown in Figures 51 and 52 for Data Sets B and C using each 
identification technique. 

The correlations achieved have standard deviations that compare favor- 
ably in level with previous APRS values, so that this technique is not 
unreasonable. The best identification technique would appear to be number 2, 
based upon Data Set C with its many combined types of patterns. This is due 
to the previously noted trend of all pattern surge airflow losses appearing 
as circumferential patterns except the pure hub-radials. Identification 
technique number 2 tends to make most combined patterns appear as circumfer- 
ential, thus providing the best correlation of AWS. 

These airflow correlations do not represent an end, but a step in the 
correlation of APRS defined at constant airflow. The method used to combine 
APRS defined at constant corrected speed (APRSN) and AWS to correlate APRS 
defined at constant airflow (APRSW) is shown in the sketch in Figure 53. 

The following items are required known values: 

• N//e^, pattern identity, distortion level 

• Clean surge pressure ratio versus corrected speed 

• Clean surge corrected airflow versus corrected speed 

• Clean surge pressure ratio versus corrected airflow 

• ^^^RSNcalC* ^^^calC 

These are all typical outputs from a compressor distortion test series. 

They permit calculation of APRSWqalC* defined in Figure 53. 
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The results of these correlations are shown in figures 54 and 55. 

Again, they represent Data Sets R and C for pattern identification techniques 
1 and 2, These correlations have standard deviations that are greater than 
those obtained with the APRSN correlations of Figures 45 and 46, Both 
pattern recognition techniques produced similar levels of standard deviation. 
This is because the APRSN correlations were best with technique number 1, 
while AWS correlations were best with technique number 2, resulting in APRSW 
correlations that are similar with both techniques. 


DIRECT CORRELATION OF APRSW 

At this point, it was desired that a correlation of loss of surge air- 
flow at constant corrected airflow be done. This was implemented with the 
same equations defined for the DIDENT approach, namely: 

P(min, 60®)^ 

^2 

However, such an approach no longer uses the parallel-compressor model 
as a basis because of the definition of APRSW. Because of this basic 
difference, the methodology for this effort will be referred to as Method E 
to discriminate it from DIDENT, Because Method E is not related to the 
parallel-compressor model, the values are now simply empirically determined 
sensitivities. The curves defining Kj[ are shown in Figures 56 and 57 for 
Data Sets B and C, respectively. There are some aspects of these curves worth 
noting. First, in Figure 56, for Data Set B, an effect of the radial extent 
of pure tip-radial patterns is seen. This was not the case with the previous 
sensitivities for APRSN. Also, the partial tip- and hub-radial patterns 
appear to fit the circumferential sensitivity curve better than they fit 
their respective pure radial sensitivity curves. The data in Figure 57, Data 
Set C, support the effect of partial hub-radial patterns appearing more like 
circumferential patterns at low corrected speeds. However, at -100% speed, 
their sensitivity is different from both pure hub-radial and circumferential 
patterns. No pure tip radial patterns exist in Data Set C so that no con- 
clusions as to extent effects can be drawn. The partial tip- radial patterns 
do not have sensitivities similar to circumferential patterns, as did Data 
Set B. Also, Kc ^ 1, shows the deviation from parallel— compressor theory. 

In general, the radial sensitivities in Data Sets B and C are widely 
different from one another for APRSW and APRSN correlations. These differ- 
ences are responsible for the apparent engine- to-engine variations previously 
noted when using the DIDENT approach. 

The pure circumferential, tip- radial (not accounting for extent) and 
hub- radial sensitivity curves in Figures 56 and 57 were used to develop a 
correlation for APRSW. They are shown in Figures 58 and 59. Again, Data 
Sets B and C with pattern recognition techniques 1 and 2 are shown. The 
standard deviations obtained for Data Set B are very nearly the same as 
those obtained in Figures 54 and 55, where APRSW was derived from APRSN and 
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AWS. This is due to the large number of pure patterns in Data Set B, In 
Data Set C, the direct correlation of APRSW is better than the derived 
correlation for pattern recognition technique number 1, and worse for 
technique number 2, 

All of these correlations can best be reviewed by referring to the 
following table of standard deviations. 


Table VI 


Comparison of Methodology Results 


Data Set B Data Set C 


Parameter Pattern 
Recognition 

No. 1 

No. 2 

Manual 

No. 1 

No. 2 

Manual 

APRSN 

0.0215 

0.0220 

0.0170 

0.0293 

0.0309 

0.0259 

AWS 

0.0160 

0.0134 

— 

0.0155 

0.0162 

— 

APRSN + AWS^APRSW 

0.0299 

0.0268 

— 

0.0360 

0.0360 

— 

APRSW, Method E 

0.0309 

0.0263 

— 

0.0344 

0.0349 

— 

APRSW, Method D 

— 

— 

0.0185 

— 

— 

0.0 360 


The values of standard deviation from the GE Method D analyses have been 
included as a reference. Interestingly, Data Set C correlates about as well 
with any one of the three approaches using APRSW. This is considered more sig- 
nificant than the standard deviations of Data Set B with its many pure 
patterns. Also, the correlations of APRSN are consistently better than the 
correlations of APRSW in the Method E analyses. Based upon the correlation 
standard deviations, the APRSN approach certainly is best. However, the 
utility of this approach to engine programs must be considered in some detail 
prior to any recommendation. 

The basic application of any distortion methodology is in the stability 
assessment, or audit, of a compression system. A methodology must be able to 
relate distortion level, sensitivity, and loss of surge pressure ratio. 

Further, loss of surge pressure ratio must be related to surge margin because 
a stability audit is the summation of all internal and external destabilizing 
influences to arrive at the amount of excess margin (or additional margin 
required) at a given operational condition. This process must be available 
to the designer during early phases of engine design so that an optimized 
target surge line, operating line, and control system may be defined. During 
this phase of the design, no design component maps may be available. 

Certainly, no test maps will be available. 

Because of the fluid nature of all parts of the engine design in its 
early stages, a APRSN value would be extremely difficult, if not impossible, 
to use in the design optimization. This difficulty precludes the recommenda- 
tion of a APRSN approach to distortion correlation. 
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The DIDENT correlations of APRSW, as seen in Table I, are rather com- 
parable in their standard deviations. Because of this, the A?RSN + AWS'> 
APRSW approach cannot be recommended due to its complexity as compared to the 
direct empirical correlation of APRSW. The recommended DIDENT formulation is 
the direct empirical correlation of APRSW, Method E, 


AP PLICATION OF METHOD E TO AN INDEPENDENT DATA SET 

The Method E concept, as used in this contract, was developed from one 
engine - the J85-13. To show that the concept is viable, it is necessary to 
successfully correlate an Independent data set from another compression 
system. This is the primary requirement for Data Set E, Data Set E is 
derived from a General Electric test of a low pressure compressor component 
and consists of 36 surge points. These surges and the various screens are 
described in Appendix I. The differences in design between the J85 compressor 
and the low-pressure compressor test vehicle are shown in Table VII. 


Table VII 


Compressor of J85— 13 and Low-Pressure Compressor 


Parameter 

J85-13 

Low-Pressure 

Compressor 

Design Speed, rpm 

16500 


13266 



Design Pressure Ratio 

~6.8 


4.1 



No. Stages 

8 


3 



Variable Stators 

IGV’s 


IGV's 

SI, S2, 

S3 

Rotor Aspect Ratios 

2.64, 

3.51, 3.58 

1.73, 

1.598, 

1.47 


3.24, 

3.1, 2.89 





2.79, 

2.27 





• SI, S2 only have inner platforms • All stators have inner 

platforms 

There are many design differences between these two machines so that a 
successful correlation of Data Set E should prove the validity of Method E. 

The screens contained in Data Set E include 180° extent one per revolu- 
tion, 50% area tip- and hub- radial, and simulated flight patterns. Other 
parameters varied in this data set include component builds and stator 
rigging. These have been included so that effects of aerodynamic variations 
within the test vehicle on Method E sensitivities can be observed. 

Data Set E does not include any variation in extent of the pure circum- 
ferential screens. A direct evaluation of ©crit possible because of 

this. The value of 6crit determined from the flight patterns to be 60°. 
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The data were then analyzed to define the curves . These are shown 
in Figure 60. These curves indicate only ininor differences in between the 
builds and stator rigging variations. This is different from the behavior 
of the J85 compressor in Data Sets B and C where significant variations in Ki 
were seen for the various engines tested. All the pure screens show consis- 
tent trends in the variation of Ki with corrected speed. The combined screens 
are less straight forward . Screen 303 indicates a sensitivity typical of a 
tip-radial pattern at 100% speed while it is typical of a circumferential 
pattern at 105% speed. Screen 402 is midway between tip- radial and circum^ 
ferential values of Ki. Both screens are combined tip- radial and circumfer- 
ential patterns . 

This behavior of the combined patterns has given further insight into 
the pattern recognition problem. Both recognition techniques previously 
used would fail to correlate these screens. NASA-LeRC has iised a different 
technique with some success on Data Sets B and C. This technique is to 
compute the following: 


APRS 


APRS 


APRS 

APRS 


T 


H 


C 


CALC 


1 - 


1 - 


P . . . 

mxn, tip rxng 

h 


min, hub ring 

Pz 




1 - 


min rake 


K. 


= Max of {APRSrj, APRSh, APRS^.} 




For a pure circumferential screen pattern, the minima are approximately equal; 
that is: 


P P P 

min, tip ring _ min, hub ring _ min rake 

?2 ?2 P2 

In this case, the APRSj^ that is maximum will be the APRSj[ associated with 
the minimum Kj^. Referring back to Figure 60, for N//6 > 90%, K^ minimum = Kq, 
for N//d < 90%, Ki minimum = Kp. This would result in circumferential 
screens having an erroneous APRSq^C < 90%. To eliminate this, the 

following logic was installed in the NASA-LeRC pattern identification tech- 
nique: 
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ir 


13 

J 


If 


P - P 

min, hiib min, tip 

h 


< 0.025. APRS^^ = APRS^ 


Otherwise APRS^^^ = Max of {APRS^, APRS^^, APRS^} 


The resulting correlation for Data Set E is shown in Figure 61. Only 
two points are significantly far from the ± 0.02 band. These are the two 
points from screen 402 that lie midway between tip-radial and circumferential 
sensitivities in Figure 2. Neither tip-radial or circumferential sensitivity 
will adequately describe these points. It should be noted that screen 402 
contains a solid plate. Past experience with screens containing a solid 
plate has been that test results are difficult to interpret and correlate 
with surge pressure ratio loss. 

The successful correlation in Figure 61 demonstrates that the Method E 
concept is a viable approach to describing surge pressure ratio loss due to 
inlet distortion for a wide range of compression system design variables, 

A very significant observation derived from this work is that the two 
separate builds and stator rigging variations have minimal effect on the 
correlation capability of Method E. This is in direct opposition to the 
J85-13 in Data Sets B and C. One must now pose the question of which data 
set is typical of all engines. If the J85-13 in Data Sets B and C is not 
representative, one of the significant potential limits to the applicability 
of Method E can be removed. 


APPLICATION O F MET H OD E T O UNSTEADY DATA 

Data Set A has been analyzed using the recommended formulation of Method 
E. The sensitivities used were those derived from the analysis of Data Set B, 
shown in Figure 56. This is consistent with the process used for the Method D 
analysis of Data Set A. 

The analysis was performed in the Stability Measurements Anlysis Lab- 
oratory (SMAL) on the analog computer. A block diagram of the analog pro- 
gram is shown in Figure 62. The analysis used two pattern recognition 
techniques. One was the NASA-LeRC technique; the second was that technique 
referred to as No. 1 in the development of the DIDENT formulation. The 
objective of the analysis was to determine a "critical time" for the J 85-13 
compressor using Method E. The data were analyzed at frequencies correspond- 
ing to 1, 1/2, 1/4, and 1/8 rotor revolution. Waveforms from data point 148 
are shown in Figure 63. 
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As in the analysis of Data Set A with Method D, peak values of APRS 
were chosen from the 20-millisecond time period prior to compressor surge. 
These values are shown in Figure 64 where APRScj^,c is plotted versus averag- 
ing time for nine of the data points . Data points 214 and 219 were excluded 
from this analysis due to their inordinately low APRS AVAIL values. 

Study of Figure 64 will show that most of the resulting APRScALC values 
are quite low compared to APRSavAIL compared to APRSqaLC by Method D. 

The results would indicate that the compressor's "critical time" is something 
less than 0.125 rotor rev. This is so far from all previous experience that 
such a conclusion is highly suspect. Furthermore, the curves in Figure 64 
for the Method E analysis are rather flat, indicating that it is highly 
improbable that the ^ERSavAIL levels could be achieved at any averaging time. 

The data were examined in detail to find a possible explanation for this 
discrepancy. The analog program was verified - the computed pressure defects 
were correct. The program logic was investigated to determine how often the 
pattern recognition technique switched from one type of pattern to another. 
The pattern choice was always stable for the 20 milliseconds prior to surge. 
For some data points, it was stable for -1 second prior to surge. These 
investigations showed that the analog analysis was mathematically correct. 


The only variable left in the formulation is sensitivity. An assumption 
was made that sensitivity is also a function of pressure defect level. To 
verify that assumption. Method D and Method E were related in a simplistic 
fashion. Consider a pure 180“ square wave circumferential distortion pattern 
so that: 


APRS 


METH D 


APRS 


METH E 



P 

1 lain 

p. 

face 


r 



SF=EF=MR=1 . 0 


Then, equating APRS^XH D j- 


. AP _ , _ min 
^ " c ^Face 
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ring av Face 


K, 
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Face 


= 1 - 
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> 

Face 


K, 
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K. 
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P . 

mxn 


Face 


1 - K. 
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min 


Face 




K, 


1 


P . 
min 


^Face 


(^~h) + '^c 


With this relationship between Method D and Method E, it can be seen that, 
for a constant value of Kq, the Method E sensitivity, Kc» would be expected 
to be an inverse function of pressure defect level. Figure 15, for example, 
shows Kc (Method D) to be constant with defect level. 

Similar relationships can be derived for radial sensitivities. The primary 
sensitivities in Method D are then relatable to Method E sensitivities. 

With the Method E sensitivities shown to be functions of pressure defect 
level, along with corrected speed and pattern, Data Sets B and C were 
reinvestigated to determine this functional relationship. For this work, 
it was assumed that Data Sets B and C could be used as a single data set. 

For each corrected speed and type of pattern, sensitivity was plotted as a 
function of pressure defect level. The results are shown in Figure 65. The 
data in this figure from Data Set A were plotted by assuming that 0.5 rotor 
revolution was the correct "critical time" so that a sensitivity could be 
computed. Study of Figure 65 shows that Data Sets B and C are not unique, 
and can be considered as one data set. Further, Data Set A pressure defect 
levels are not all contained within the steady-state data. In order to 
analyze Data Set A using this DIDENT approach, one would need to extrapolate 
the sensitivity curves to defect levels typical of Data Set A. The results 
would be a function of the extrapolation chosen. For this reason. Data Set 
A was not further analyzed. 

The most important observation is that Data Sets B and C can be con- 
sidered as one data set. This explains the apparent engine- to- engine varia- 
tion of Method E sensitivities - it is just a result of higher distortion 
levels in Data Set C compared to Data Set B. One of the limits to the use 
of Method E is then removed. However, the number of patterns required to 
define sensitivites has been increased because of the need for a range of 
pressure defect levels. Further, the distortion levels used in defining 
sensitivities must encompass those anticipated in the engine installation. 

It is now possible to correlate Data Sets B and C as one. This has 
been done in Figure 66 using both pattern recognition techniques developed 
herein. This will allow a comparison of this correlation with the previous 
individual correlations of Data Sets B and C. Again, the comparison is best 
made with the standard deviations. These are shown in Table VIII where 
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2 B+C is the summed standard deviations from individual correlations using 
the number of data points in each data set as weighting factors . 


Table VIII 


Comparison of Data Sets B and 



Pattern 

Recognition 

Data Set 

No. 1 

No. 2 

D 

B 

0.0309 

0.0263 

0.018 

C 

0.0344 

0.0398 

0.036 

E B+C 

0.0324 

0.0332 

0.027 

B and C 

0.0295 

0.0,309 

— 


The standard deviation of the combined data sets shows a better corre- 
lation than the summed standard deviation of each data set for both pattern 
recognition techniques. In the case of pattern recognition technique 
number 1, the combined data sets have a smaller standard deviation than 
either of the individual data sets. 

Method E can be described as a viable distortion methodology that can 
be applied to any compression system. The sensitivities in the recommended 
formulation used here are functions of speed, pattern, and pressure defect 
level. They are not functions of an individual engine. The number of 
patterns needed at each speed to define the Method E sensitivities is esti- 
mated in Table IX, 


Table IX 

Estimated Patterns Needed to Implement Method E 
Pattern Pressure Defect Levels Sensitivity 


180 one-per-rev Low, 

120° one-per-rev Low 

90° one-per-rev Low, 

60° one-per-rev Low 

45° one-per-rev Low 

Hub-radial, 40% extent Low, 

Hiib- radial, 20% extent Low, 

Tip-radial, 40% extent Low, 

Tip- radial, 20% extent Low, 


high 



high 

I 

/ crit ^ 

medium, high 
high 

J 

) Kh 

medium, high 
high 


) % 


A total of 17 patterns at each speed would be needed to implement 
Method E. The many radial screens are included so that effects of extent 
of the radial defect may be observed. Note that one screen can typically 
produce acceptable patterns over a range of speeds. This does not say 
that 17 screens are needed at each speed. 
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CONCLUSIONS 


The development of the General Electric distortion methodology. Method D, 
has been documented. Method D has been applied to steady-state data to show 
its capability to correlate inlet pressure distortion data to loss of surge 
pressure ratio. The sensitivities used in Method D that were derived from 
the steady-state data were used to analyze unsteady data to determine the 
"critical time" of the J85-13 compressor. The "critical time" is 0,5 rotor 
revolution, as determined by these data. 

It was shown, using Method D equations, that output averaging will produce 
pressure defects greater than or equal to input averaging. However, due to the 
nonlinear sensitivities in Method D, it can only be stated that output aver- 
aged values of surge pressure ratio loss are not equal to input averaged values 
of surge pressure ratio loss. The use of output averaged pressure defects 
to compute surge pressure ratio loss can then lead to highly erroneous values , 

Analog and digital analyses of unsteady data have been compared. Analog 
analysis was shown to agree with digital analysis within jjO.Ol of the Method D 
pressure defect terms. Within the ground rules established for this comparison, 
analog analysis is faster and less costly than digital analysis. 

Two formulations of the NASA-LeRC distortion methodology, DIDENT, were 
investigated. One formulation correlated the loss of surge pressure ratio 
at constant corrected speed. The other correlated the loss of surge air- 
flow at constant corrected speed and combined this with the previous correla- 
tion to obtain a correlation of the loss of surge pressure ratio at constant 
corrected airflow, A methodology that used DIDENT equations to directly 
correlate the loss of surge pressure ratio at constant corrected airflow was 
implemented. It was concluded that this last formulation. Method E, provided 
the most useful results. 

The use of DIDENT and Method E requires that some form of pattern recog- 
nition be employed to discriminate between hub-radial, tip-radial, and circum- 
ferential patterns. Two techniques were developed. Both were based on 
pressure defect levels and locations in the entrance plane. Resulting corre- 
lations showed that neither technique was completely satisfactory. Work at 
NASA-LeRC on this problem in their work with DIDENT showed that a technique 
that maximized the computed loss of surge pressure ratio was a more accurate 
technique. 

Method E was applied to an Independent data set to verify that the 
formulation was applicable to data other than that set used to develop it. 

The pattern recognition technique used was the NASA-LeRC process. The results 
showed a good correlation of loss of surge pressure ratio. This would indi- 
cate that Method E is a viable distortion methodology, 

A set of unsteady data were analyzed using analog techniques with the 
Method E formulation. This analysis was intended to determine the "critical 
time" of the J85-13. The analysis, instead, showed that the Method E sensi- 
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tivities are functions of pressure defect level. Because the unsteady data 
defect levels were typically higher than steady-state data, it was decided not 
to analyze the unsteady data for "critical time." Such an analysis would 
require an arbitrary extrapolation of sensitivities to higher defect levels. 
This effect, however, was used to explain the apparent engine-to-engine vari- 
ation seen in Method E sensitivities - it is just an effect of pressure defect 
magnitude. 

Method E can then be summarized ; 

• The Method E formulation is derived from the DIDENT methodology 
developed by NASA-LeRC, DIDENT utilizes the parallel-compressor 
model as its basis. Method E, however, is considered to be an 
empirical formulation not restricted by any parallel compressor 
model assumptions. 

• The correlation capability of Method E was shown to be comparable 
to Method D, Also, the apparent engine-to-engine variation in the 
Method E sensitivities has been shown to be nonexistent. 

• At this point in its development. Method E requires that the follow- 
ing parameters be defined to implement the methodology: 

• ®crlt “ angle over which the minimum pressure is integrated. 

• Kq ” circumferential sensitivity, a function of corrected speed 
and distortion level. 

• Kjj - hub-radial sensitivity, a function of corrected speed, dis- 
tortion level, and (potentially) of extent of low pressure 
region. 

• Kx ” Tip-radial sensitivity, a function of corrected speed, dis- 
tortion level, and (potentially) of extent of low pressure 
region. 

The sensitivities would need an estimated 17 patterns at each speed to be 
defined. Again, it is noted that one screen can produce acceptable patterns 
over a range of corrected speeds . 
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APPENDIX I 
DATA SET E 


Data Set E consists of 36 surge points whose speeds and APRS^y^jL values 
are tabulated in Table I-I, Instrumentation used to define the patterns is 
shown in Figure 67. Screen patterns are shown in Figures 68 and 69, 
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Table I-I 


Data Set E 





Stator 

APRS 

APRS 

SCN 

SPD 

Build 

Schedule 

AVAIL 

CALC. 

100 

105 

1 

Nom, 

0.0941 

0.0771 


100 

1 

Nom, 

0.0681 

0.0615 


95 

1 

Nom, 

0.0390 

0.0349 

101 

95 

1 

Nom, 

0.0969 

0.0988 


90 

1 

Nom, 

0.0480 

0.0536 


85 

1 

Nom, 

0.0317 

0.0307 


75 

1 

Nom, 

0.0 

0.0006 

201 

105 

1 

Nom, 

0.0360 

0.0360 


100 

1 

Nom, 

0.0394 

0.0413 


95 

1 

Nom, 

0.0575 

0.0548 


90 

1 

Nom, 

0.1029 

0.0842 

202 

95 

1 

Nom, 

0.0965 

0.0975 


85 

1 

Nom, 

0.0822 

0,0822 


75 

1 

Nom. 

-0.0094 

-0.0094 

206 

105 

1 

Nom. 

0.0371 

0.0371 


100 

1 

Nom. 

0.0363 

0.0363 


95 

1 

Nom, 

-0,0517 

-0.0517 

401 

105 

1 

Nom, 

0.0060 

0.0098 


100 

1 

Nom, 

0.0 

0.0147 


97.5 

1 

Nom. 

0.0024 

0.0137 

402 

90 

1 

Nom, 

0.0805 

0.1225 


85 

1 

Nom. 

0.0632 

0.0911 


75 

1 

Nom, 

0.0093 

-0.0036 

100 

105 

4 

Nom, 

0.0667 

0.0592 


105 

4 

+5" 

0.0711 

0.0571 


100 

4 

Nom, 

0.0689 

0.0552 


100 

4 

-5*=* 

0.0739 

0.0542 


95 

4 

Nom, 

0.0453 

0.0361 

201 

100 

4 

Nom, 

0.0279 

0.0269 


100 

4 

-5*’ 

0.0256 

0.0270 


95 

4 

Nom, 

0.0453 

0.0470 

303 

105 

4 

Nom, 

0.0485 

0.0235 


100 

4 

Nom, 

0.0262 

0.0262 


100 

4 

Nom, 

0.0221 

0.0272 


95 

4 

-7° 

0.0263 

0.033 


95 

4 

Nom, 

0.0263 

0,0294 
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LIST OF SYMBOLS 


Method A 


e 


0 ) 

n 

•^n 

a 


Pf 

q 

Ke 

B 

b 

Kr 

aprs^est 


“ total pressures in a ring 

- Fourier coefficients 

- angular position around inlet 

- Fourier coefficient weighting function, empirically determined 

- nth Fourier coefficient 

- amplitude coefficient 

- diameter of jth ring 

- circumferential ring diameter weighting function, empirically 
determined 

-• face average pressure 
-jth ring average pressure 

- dynamic head 

- circumferential index 

- radial ring diameter weighting function, empirically 
determined 

- radial speed function, empirically determined 

- radial index 

- loss of surge pressure ratio measured from engine/component 
tests 
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Method B 



- ith pressure in a ring 

- jth ring average pressure 

- Fourier series coefficients 

- angular position around face 

- sum of harmonic amplitudes for jth ring 

- nth harmonic of Fourier series 

- multi-per-rev weighting function, empirically determined 

- empricial ring weighting function for circumferential 
distortion 


- circumferential index, empirical 

- circumferential sensitivity, empirical 

APRS^eST " loss of surge pressure ratio measured from engine/component 
test 


- loss of surge pressure ratio due to circumferential distortion 

Pp “ face average pressure 

Wj. - empirical ring weighting function for radial distortion 

- sum of harmonic amplitudes of radial distortion 

- radial distortion sensitivity, empirical 

- loss of surge pressure ratio due to radial distortion 

NDrji - loss of surge pressure ratio for generalized pattern 
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Method C 


P, 

1 

P ■ “ 

■^min 

Pf 

n - 

£0 - 
IDC 

Kc 

APRS cir cum” 

Ki 

k 

IDR-x 

IDRh 

Kr^ 

APRStoTAL " 
APRS test “ 


average pressure in ith ring 
minimum pressure in ith ring 
face average pressure 

number of low pressure regions around each ring 
empirical weighting factor on n 
circumferential index 
circximf erential sensitivity 

loss of surge pressure ratio due to circumferential distortion 

empirical ring interactor factor 

radial extent factor 

tip radial index 

hub radial index 

tip sensitivity 

hub sensitivity 

loss of surge pressure ratio 

loss of surge pressure ratio 
pattern 

loss of surge pressure ratio 
tests 


due to radial distortion 
due to general distortion 

measured from engine /component 
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Method D parameters are defined in Figure 5 


DIDENT AND Method E 

P3/P2)d 

?3/P2)c 
(Pmin. 60°)2,r 

K± 


- distorted compressor pressure ratio 

- undistorted compressor pressure ratio 

- minimum pressure over 60° sector in plane 2,D (inlet 
plane), hub ring or tip ring only 

- sensitivity, empirically determined. Related to 
parallel-compressor pressure ratios for DIDENT. Con- 
sidered to be an empirical coefficient only for 
Method E 
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FIGURE 1 - DEFINirrON OF LOSS OF SURGE PRESSURE RATIO AND LOSS OF SURGE AIRFLOW 
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PROCEDURE FOR CIRCUMFERENTIAL CONTRIBUTION 
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PROCEDURE FOR RADIAL CONTRIBUTION 
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FIGURE 2 - METHOD A PROCEDURK 



PROCEDURE FOR CIRCUMFERENTIAL CONTRIBUTION 
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PROCEDURE FOR RADIAL CONTRIBUTION 
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FOURIER SERIES APPLIED TO AVERAGED RING PRESSURES 
P = P, 
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FIGURE 3 - METHOD B PROCEDURE 
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iPRS 

APRS TEST 

TEST 


PROCEDURE FOR CIRCUMFERENTIAL CONTRIBUTION 
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PROCEDURE FOR RADIAL CONTRIBUTION 




• FOR EACH RING DETERMINE LEVEL FACTORS 



• DEFINE TIP AND HUB LEVEL FACTORS 

TDRT ? i IDR. . + IDR, , 

TDRH = Jl IDR, + IDR . „ 

1=1 1=2 

= f (rndial extent, N/VQ) 

• FOR BASIC RADIAL APRS^^^^ DETERMINE SENSITIVITY 



IDR„ 


'•‘^ 1,5 

DEFINE RADIAL INDEX 


PROCEDURE FOR COMDINEF) PATTERNS 

• A FR^ TOTAL = ^P^CIRCUM + 


Aprs rad 


FIGURE 4 - METHOD C PROCEDURE 
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METHOD D CIRCUMFERENTIAL DISTORTION 


FIGURE 5 - METHOD D PROCEDURE 
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FIGURE 8 - METHOD A - GE-FIOO BU/7 & IGV'S & O.S, 




FIGURE 9 - METHOD A - GE-FlOO BU/7 & IGV^S FIGURE 10 - METHOD A - GE-FlOO BU/6 
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FIGURE 11 - METHOD A NASA - LeRC DATA SET B 



FIGUBIE 12 - METHOD B CORRELATION 
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BOUNDS OF NASA - J85 DATA- 

cr=.0234 I I 
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FIGURE 16 - DATA SET B CIRCUMFERENTIAL FIGURE 17 - DATA SET B EXTENT FUNC- 
SENSITIVITY VS ANGULAR TION VS ANGULAR EXTENT 

EXTENT 



N~PER REVS 


FIGURE 18 - MULTI-PER REV FUNCTION VS NUMBER OF LOW PRESSURE ZONES 
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FIGURE 23 - J -85 DATA SET 


B SUPERPOSITION FACTOR VS. ^PRSj^/ApRSq 




FIGURE 24 - J-85 DATA SET B CORRELATION 
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FIGURE 25 - CIRCUMFERENTIAL SENSITIVITY AS A 
FUNCTION OF CIRCUMFERENTIAL DIS- 
TORTION FOR VARIOUS ANGULAR EXTENTS 
(METHOD D+) 
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FIGURE 26 - CIRCUMFERENTIAL EXTENT FUNCTION 
VERSUS ANGULAR EXTENT (METHOD 
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FIGURE 27 - CIRCUMFERENTIAL SENSITIVITY AS A 
FUNCTION OF 180'* - NORMALIZED 
CIRCUMFERENTIAL DISTORTION FOR 
VARIOUS RADIAL EXTENTS (METHOD D+) 
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FIGURE 28 - RADIAL EXTENT FACTOR VERSUS RADIAL 
EXTENT (METHOD D+) 
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FIGURE 29 - RADIAL SENSITIVITY AS A FUNCTION 
OF RADIAL DISTORTION (METHOD D+) 




- SUPERPOSITION FUNCTION VERSUS RADIAL 
LOSS IN SURGE PRESSURE RATIO 
(METHOD D+) 
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FIGURE 31 - MEIHOD D+ CORRELATION, DATA SETS B AND C. 
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FIGURE 33 - SURGE NO. 141, DATA SET A 
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AVERAGING TIME, ROTOR REVOLUTIONS 


FIGURE 34 - HISTOGRAM OF AVERAGING TIME, 
METHOD D, DATA SET A 
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FIGURE 35 - STANDARD DEVIATION VERSUS AVERAGING 
TIME, METHOD D, DATA SET A 
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FIGURE 38 - DATA SET A ANALOG ANALYSIS BLOCK DIAGRAM 
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FIGURE 39 - data set A, POINT 82 
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FIGURE 4o - COMPARISON OF DIGITAL AND ANALOG ANALYSES DATA SET A 
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FIGURE 45 - CONSTANT CORRECTED SPEED CORRELATION OF DATA SETS B AND 
C USING PATTERN RECOGNITION TECHNIQUE NO. 1 
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FIGURE 46 - CONSTANT CORRECTED SPEED CORRELATION OF DATA SETS B AND 
C USING PATTERN RECOGNITION TECHNIQUE NO. 2 
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• PARALLEL - COMPRESSOR SECTORS 
■ PREDICTtD SURGE POINT 


FIGURE 47 - PREDICTION OF SURGE AIRFLOW USING A 
PARALLEL - COMPRESSOR MODEL 
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FIGURE 53 - SUMMATION OF ZEPRSN AND AWS INTO APRStf 
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FIGURE 55 - C0RRELA.TI0N OF ApRSN AND AWS INTO ApRSW, PATTERN 

RECOGNITION TECHNIQUE NO. 2 
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FIGURE 56 - METHOD E DENSITIVITIES FOR LOSS OF 
SURGE PRESSURE RATIO AT CONSTANT 
AIRFLOW CORRELATION, DATA SET B 
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FIGURE 57 - METHOD E SENSITIVITIES FOR LOSS OF 
SURGE PRESSURE RATIO AT CONSTANT 
AIRFLOW CORRELATION, DATA SET C 
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FIGURE 60 - DIDENT SENSITIVITIES AT CONSTANT CORREC 
LOW PRESSURE COMPRESSOR RIG TEST - DATA £ 
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FIGURE 6l - DIDENT CORRELATION OF SURGE PRESSURE RATIO LOSS AT CONSTANT CORRECTED 
AIRFLOW FOR LOW PRESSURE COMPRESSOR RIG TEST DATA SET E 
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FIGURE 65 - METHOD E ANALYSIS WAVEFORMS, POINT l48, DATA SET A 
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FIGURE 64 - CALCULATED PEAK APRSW AS A FUNCTION OF AVERAGING TIME 
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FIGURE 66 - CORRELATION OF DATA SETS B AND C AS ONE DATA SET 
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